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ABSTRACT 

Traditional analysis oi g = 2(1 + ^) experiments for charged leptons use a classical spin 
vector picture. For muons, we here employ a more exact Dirac quantum four component 
spinor theory. Survival probabilities (including wave packet effects) are computed. These 
oscillate with the frequency O = (KeB/Mc) as has been assumed in previous muon {g — 2) 
experimental analyses; i.e. muon survival probability oscillations are already at the root 
of previous succesful {g — 2) measurements. Further oscillations should also be observed 
if mixed neutrino mass matrices were to enter into the reactions (e.g. tt"*" —>■ ii'^ + v^j) 
producing muons. 
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1. Introduction 



Measurements of the charged lepton anomalous magnetic moments [1-5], i.e. the g- 
factor 

^ = 2(l + «), (la) 

are of considerable interest. The calculation of k as a function of the coupling strength 
a = (e^/hc), 

gives one some confidence (when compared with experiment) that perturbative quantum 
electrodynamics (and perhaps some other gauge theories) makes sense. New precision 
measurements of k for the muon have been proposed[6] for obtaining additional insights 
into strong and weak inter actions [7]. 

The analysis of {g — 2) measurements often proceeds from classical equations [8, 9]. For 
example, a charged lepton is considered to follow a classical orbit in a magnetic field 

Even the spin- 1/2 pseudo- vector is considered to follow a classical equation of motion 



dr ~^ Mc ' ^"^^ Mc3 ' ^ 



Such classical spin equations can serve only as a rough guide to the inherent quantum spin 
interference exhibited in laboratory high precision measurements of {g — 2). 

Consider the precision measurements of (5^ — 2) for the muon [2-5]. Muons (produced 
by the decay (say) 7r~ — > iJ,~ + 9^) are injected into a ring with a uniform applied magnetic 
field B. The experimental quantum survival probability of the muon in the ring to decay, 
(via fx~ ^ e~ + Ue + J-'^ with a detected electron energy above a threshold value) has been 
fit to the theoretical functional form 

p u\ _ „-Mc^rt/£ ^ 1 + ^ cos(Ot + cfy ^ 

/^^-^e-+.,+.eW - e i 1 + Acoscj> >■ 

In Eq.(3), M, and represent (respectively) the mass, energy, and intrinsic lifetime 
of the muon; t is time in the laboratory reference frame, and 

n = (^) (4) 

is the experimental frequency measured in the muon survival probability due to quantum 
mechanical amplitude interference. The frequency in Eq.(4) is central for the experimental 
determination of k in Eq.(l). Eq.(3) for the muon survival probability (as determined by 
a detected energetic electron) is central for the {g — 2) measurement. Eq.(3) represents a 
quantum mechanical charged lepton oscillation which follows from the Dirac equation and 
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cannot be fully understood in purely classical terms. This point will be discussed in detail 
in the work which follows. 

Starting from the work of Schwinger[10,ll], the magnetic moment anomaly parameter 
K was defined by the manner in which the vacuum renormalized lepton mass depended on 
an applied magnetic field. For example, the Dirac-Schwinger equation for a muon moving 
in a magnetic field has the form 

( - in-f^'d/^ + Mc- i{nr/2c))i;{x) = 0, (5a) 

where the gauge derivative is defined as 

eA 

d^, = df,-i{^), F^, = d^A, - d^A^, (56) 
and the renormalized mass (matrix) has the form 

Our purpose is to discuss in detail the notion of charged lepton oscillations as they 
have appeared in previous {g — 2) measurements employing the quantum mechanical Dirac- 
Schwinger Eqs.(5), rather than the classical Eqs.(2). The quantum mechanical superposi- 
tion of amplitudes viewpoint is by far the more fundamental. In Sec.2, the Dirac-Schwinger 
Eqs.(5) will be derived from the usual charged current definition of k, i.e. 

J^(a;) = eciPixh^'iPix) + {^)d.{iP{x)a^''il;{x)). (6) 

It will be shown that the second term on the right hand side of Eq.(6) gives rise to the 
second term on the right hand side of Eq.(5c). The exact energy spectrum of a charged 
lepton (with k ^ 0) in a uniform magnetic field will be derived in Sec. 3 from the Dirac 
equation. In Sec. 4, the survival amplitude for a muon in a magnetic field will be computed 
from solutions of the Dirac equation and the experimental modulation frequency Eqs. (3) 
and (4) will be derived. The general nature of charged muon oscillations (including those 
induced by neutrino oscillations) is discussed in the concluding Sec. 5. 

2. The Dirac-Schwinger Equation 

The action for a spin- 1/2 charged particle may be written as 

S = J d^xi^iih-^i'd^, - Mc)^ + ^ / d'^^J^^i^- (7) 
From Eqs. (6) and (7) it follows that 

S = j d^xi;{inrd^ - Mc)i; - {^^) J d^xi^a^^^^jd^A^, (8) 
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where is defined in Eq.(5b) and integration by parts has been performed on the second 
term in Eq.(8). From Eqs.(5c) and (8), one may write the action 

from which the Dirac-Schwingcr Eq.(5a) follows. However, implicit in the above consider- 
ations is the complex mass replacement rule 

Mc^ ^ Mc^ -i{nT/2), (10) 

which will be used throughout this work to describe muon decay rates. 

For a particle in a uniform magnetic field, the mass matrix in Eq.(5c) has two eigen- 
values, 

Mw± = M±w±, (11) 
where w± are four component spinors and 

M^c^ = Mc'^{^)=Mc^^i^). (12) 

The mass splitting is thus determined by Eq.(4) as 

AMc^ = \M+ - M-\c^ = nn. (13) 

In analogy to X-meson[12] and S-meson physics, there may be a temptation to employ 
a phase shift 9 = —{c^ /h){M.+ — M.-)t = fir in terms of the proper time of the classical 
Eqs.(2). Such temptation is ill advised since the experimental phase interference in Eq.(3) 
involves laboratory time t with the phase fit and not the "proper time" r — {Mc^/E)t 
with phase Or. This kind of error does not arise if one merely solves the Dirac-Schwinger 
Eqs.(5). And this we shall now proceed to do. 

3. Energy Wave Functions 

The Dirac-Schwinger Hamiltonian in the laboratory frame in which there is an applied 
uniform magnetic field is given by 

H = ca-U + l3{Mc^ -{Kne/2Mc)a -B), (14) 

where 

n = p - (e/c) A, [n^, Uj] = i(ne/c)eijkBk, (15) 

and 

a = 75(7. (16) 

We seek solutions to 

H^{r) = EiP{r). (17) 
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With the matrix p2 defined as 

P2 = iPl5, pI = 1, (18a) 

one notes that 

Hp2 + P2H = 0. (186) 

Given a particle solution in Eq.(17) with E > 0, one may also construct anti-particle 
(negative energy) solutions 

Hp2ilj{r) = -Ep2il;{r). (19) 

Here we consider only positive energy (E > 0) eigenstates. 

Note that orbital angular momentum and spin angular momentiim about the 2;-axis 
(chosen parallel to the magnetic field) are not separately conserved; i.e. with L = r x p 
and S = ha/2, 

[H, S,] = -[H, L,] = 2zc75(n X S),. (20) 
On the other hand the total angular momentum along the magnetic field axis is conserved 

J, = L, + S„ [H,J,]=0. (21) 

Another conserved quantity is II^, 

[H,U,] = 0. (22) 

For motions in a plane perpendicular to the magnetic field Ilzip{r) — 0, i.e. there is zero 
motion along the magnetic field axis. 

Note that for the classical Eqs.(2) of motion in the plane (zero motion along the 
magnetic field axis) we have {dSz/dt) classical = 0. However, in the quantum mechanical 
treatment (dSz/dt) quantum = {i/'h)[H, Sz] 7^ 0. The classical equations of motion thereby 
contain integrals of motion not present in the physical quantum mechanical treatment. 
The point is that a particle with total spin s = {fi/2) is not a classical spinning particle. 
A classical spinning particle has a total spin angular momentum large on the scale of Ti. 

If one employs cylinder coordinates r = (p, 0, z), and the Landau gauge for the vector 
potential in B = V x A, 

Ap = Q, A^={Bp/2), ^, = 0, (23) 
then for zero motion along the magnetic field axis 

H^t/; = -in{dil}/dz) = 0, (24) 

so that 

= V'(P, 0) = H{p)e'^^ + «T(p)e^^'"-'^^- (25) 

In Eq.(25), W]^,t(p) are four component spinors and conservation of total angular momen- 
tum has been invoked; i.e. [H, J^] = allows us to choose the eigenfunctions in Eq.(17) to 
also be eigenfunctions of — + S^- Thus, the choice 

JMp, cf>) = h{m- {l/2))^|;{p, 0), (m = 0, 1, 2, ...) , (26) 

5 



determines the coherent superposition of spin up and spin down spinors in Eq.(25) which 
enter into high energy experimental measurements of (g — 2). 

From Eqs.(14), (15), (23), and (25) we find the exact energy spectrum 



^J Mlc^ + 2necBN ^ {KheB /Mbc)^ Mlc^ + 2necBN, (27a) 



where N = Up + \m\ and Up is an integer counting nodes in the radial wave functions. If 
there were a momentum component pz parallel to the magnetic field, then Eq.(27) would 
be modified to 



E±,n{Pz) = Y M|c4 + c2p2 + 2necBN ^ (nheB/MBc) Mlc^ + 2necBN. (28) 

Eq.(28) describes rigorously the exact energy spectrum of a Dirac spin-1/2 particle with a 
^-factor defined as (7 = 2(1 + k) moving in a uniform magnetic field B = |B|. 

For magnetic fields of laboratory size (say less than a few Tesla), the inequality 
nheB « M'^c^ holds by an overwhemingly large margin. Thus, very accurately Mb 
may be replaced by M, the liftime of the energy eigenstate to weak muon decays may be 
determined by 

^±,N - (^)r, (29) 

and the Bohr transition frequency 

E^^M - E+,N = nn, (30) 

where Cl is given by Eq.(4) independent of N, i.e. independent of the energy E±^n itself. 
In deriving the transition rate T±^n for the weak decay of the muon energy eigenstate, 
Eqs.(lO) and (27a) were employed with << Mc^ also satisfied again by a very wide 
margin. Within the same high accuracy approximation one may write 

E±,N+n - E±,N = {-^ ) , 1 < |n| « N. (31) 

Eqs.(29) and (30) constitute the quantum mechanical basis for Eq.(3) which lies at 
the root of the experimental analysis of the muon {g — 2) measurement. More general 
statements can be deduced about lepton polarization in magnetic fields on the basis of the 
fundamental differences between the quantum and the classical viewpoints. According to 
the classical Eqs.(2) of motion, it is quite possible for a high energy particle with E » Mc^ 
and 7^ moving in a circular orbit in a plane normal to B to have the spin polarized 
parallel to the magnetic field. For a quantum Dirac particle in a high energy (on the mass 
scale) eigenstate and with k ^ 0, the spin polarization is required by the Dirac- Schwinger 



6 



Hamiltonian to be (almost) perpendicular to the magnetic field. We have discussed above 
the reasons for this difference. For the classical equations of motion is conserved. For 
the quantum equations of motion, is not conserved. This has an effect on lepton beam 
polarizations in high energy machines, e.g. LEP. We note (in this regard) that (at LEP) 
the electron cyclotron frequency Qc = (ecB/E) is certainly not large when compared with 
the anomaly frequency Q = {kbB /mc) because even though k << 1, the energy E » mc^ 
is quite high. 

4. The Muon Survival Amplitude 

Survival amplitudes are conventionally defined as 

S{t)^<,p{0)\,p{t)>, (32) 

i.e. the amplitude that a quantum object with a wave function '0(0) at time zero will still 
be in the same state after a time t. We take the scalar product for the muon to be 

S{t) = J d^ri;\r,0)4;{r,t). (33) 
For the muon moving in a uniform magnetic field 

V'(r, ^) = XlXl (ixp{-Ts^Nt/2)exp{-iEs^Nt/n)cs^Ni^s^N{r), (34) 

s=± N 

so that 



'^W = XI XI \cs,Nfexp{-Ts,Nt/2)exp{-iEs,Nt/h). (35) 

s=± N 

The Survival probability P{t) — \S{t)\'^ is then 



Pit) = X \(^s,N\'^\^s',N'\'^exp{-Js,N,s',N't/2)cOs{u!s,N,s',wN't), (36a) 
s,N,s'N' 

where 

ls,N,s',N' = ^s,N + ^s',N', ^s,N,s',N' = -£'s,iV — Es',N'- (366) 

From Eqs.(29), (30), (31), and (36) one finds 

P{t) = ^J dW{E) J dW{E')exp{—^E,E't/'^){cos{(jjE,E't) + cos{LiJE,E't + (37a) 
where dW{E) is the probability that the muon has an energy in the interval dE, 

1E..,=T{— + —), .,,„ = (_)(_-_), (376) 
and 0, is defined in Eq.(4). Eq.(37a) is more simply written as 

P{t) = lK{t){l + cos{Qt)), (38a) 



where 



K{t) = J dW{E) J dW{E')exp{-jE,E't/2)cos{ujE,E't) (386) 

describes the spreading of the muon wave packet moving around in the magnetic field. 
Such wave packet spreading depends on the precise nature of muon energy distribution 
dW{E)/dE. For a Gaussian muon energy distribution, with mean energy E =< E > and 
deviations from the mean SE = V < E'^ > - < E « E, 

K{t) = exp{ g )exp{—^), 

Comparing Eqs.(3) and (4), which have been previously employed [2-5] for the analysis 
of muon ((7 — 2), to our central Eqs.(38) we note some similarities and some differences: (i) 
Our cos{0,t) does not have an arbitrary phase (f) as in cos{Qt + 0). The reason for this is 
that as theorists we may declare that the initial wave function is known at exactly "time 
zero" . The experimentalists actually have to build the clock and provide the beam (quite 
different from a mere theoretical pronouncement), and "time zero" has fluctuated a bit 
in past experiments. This has been reported as a fluctuating (from day to day) phase cp. 
(ii) Eq.(3) has a coefficient A < 1 while Eqs.(38) correspond to A = 1. The reason for 
this is firstly that we have assumed polarization perfectly in the plane perpendicular to 
the magnetic field and nothing in an experiment is perfect. Secondly, we have calculated 
in Eq.(38) the total survival probability while Eq.(3) refers to a muon decay ejecting 
an electron above a high energy threshold. The value of A depends at least in part on 
experimental cuts in the reported electron counts, (iii) Eq.(3) contains only the intrinsic 
exponential decay, while our Eqs.(38) also contains wave packet spreading effects modeled 
assuming a Gaussian energy distribution. The reason is that wave packet spreading effects 
have to be treated using a quantum mechanical wave (Dirac-Schwinger) equation. The 
spreading of the wave packet is in a time ts ~ (l/Suj). This should not adversely effect the 
accuracy of measured k values for SE/E < 10~^ in laboratory magnetic fields of order 1.5 
Tesla. 

5. Conclusions 

From the viewpoint of quantum mechanics, the reason for the oscillations of any total 
survival amplitude 

5tot(t)=<^'|e-^^*/'^|^'>, (39) 
is in the nature of the energy probability distribution 

dWtotiE) =< ^\S(E - H)\^ > dE, (40) 



I.e. 



Stot(t) = J dWtot{E)e-'^'/\ (41) 
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For example, a two peaked Gaussian distribution in energy (A << E), 



dWtn,o peak{E) = (l/2)Vl/27rA2 J2 exp{ - {E - E + snu;o/2f /2A^)dE, (42a) 

s=±l 

gives rise to an oscillating in time survival amplitude 

^two peak 

(t) = cos{ijJot/2)exp{-iEt/h)exp{-A^t^ /2h^). (426) 

For another example, if the total energy distribution is a convolution of two energy distri- 
butions 

dWtot{E) _ r^ dWijE + e) dW2{e) 

dE J ^ dE de ' ^ 

then the survival amplitude is the product of two survival amplitudes 

Stot{t) = Si{t)S2{t). (436) 

Thus Eq.(38a) represents a convolution of a normal weak muon decay with a two peaked 
(at E ± {hQ/2)) distribution separated by f2 = (neB/Mc), i.e. the observed muon {g — 2) 
oscillation frequency. 

Consider the notion of neutrino mass mixing[13]. In two previous experiments[14,15], 
attempts were made to find neutrino mass matrix mixing by observing a muon propagating 
in a magnetic field. One tried to measure a double peak in the muon energy distribution 
using a time scale very short compared with The point was that the neutrino mass 

from 7v~^ fj.'^ + would induce in the muon energy distribution a double peak from (say) 
two possible masses for the neutrino. These two possible masses would yield two possible 
muon energies, and thus exhibit a double peak in the muon energy distribution. Such a 
double peak induced by neutrino masses did not actually appear within the accuracy of 
the previous experiments. However, in those experiments one attempted to measure the 
muon energy distribution directly. We are here suggesting that a more sensitive probe of 
energy splitting is of the type already successful in muon {g — 2) experiments; i.e. one 
should look for the double peak in energy via the muon survival probability oscillation in 
time. The resulting extra oscillation should appear in a modified form of Eq.(3) 

P,-^e-+.^+.St) - e-^^^"/^(l - sin\29)siriHn't/2))C-±^^^), (Ua) 
where 9 is the two flavor neutrino mixing angle, and 

nn'^E^^)'. (446) 

We have assumed that the pion has decayed into a forward moving muon (in the laboratory 
frame) and neutrino with E^^ » M^^c^. The neutrino masses (differing by Am) are 
assumed very small on the scale of pion and muon masses. For a freely moving muon (i.e. 
= 0), only the neutrino induced oscillations survive as previously discussed [16]. 
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Finally, there are certainly two characteristic frequencies associated with (g — 2) 
measurements, i.e. the anomaly frequency O = (KeB/Mc) and the cyclotron frequency 
Oc = (ecB/E). The cyclotron frequency Oc has been measured directly from emitted 
radiation. To some extent, the anomaly frequency Q should also enter into the radiation 
spectrum, although the authors are not aware of any direct measurements of this sort. If 
there is also neutrino mass mixing, then a third frequency Q' enters into the problem as 
well. 
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